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Abstract 

The Higgs boson spectrum of the Next-to-Minimal Supersymmetric Standard Model is ex- 
amined. The model includes a singlet Higgs field S in addition to the two Higgs doublets 
of the minimal extension. 'Natural' values of the parameters of the model are motivated by 
their renormalization group running and the vacuum stability. The qualitative features of the 
Higgs boson masses are dependent on how strongly the Peccei-Quinn U{1) symmetry of the 
model is broken, measured by the self-coupling of the singlet field in the superpotential. We 
explore the Higgs boson masses and their couplings to gauge bosons for various representative 
scenarios. 



1 Introduction 



Supersymmetric models ^ [2] take an important step toward a solution of the hierarchy prob- 
lem [3j by stabilizing the ratio of the electroweak and Plank/GUT scales. In the Standard Model 
(SM), the Higgs boson mass gains radiative corrections which depend quadratically on the cut-off 
scale of the theory (usually taken to be the GUT scale Mqut — 2 x 10^^ GeV), threatening to 
generate a mass which is far too large to explain electroweak symmetry breaking. In order to 
stabilize this mass at a phenomenologically acceptable scale (below a TeV), one must fine tune 
the parameters of the model, creating an accidental cancellation between the bare mass and the 
quantum corrections. By introducing a new symmetry, a supersymmetry between bosons and 
fermions, one introduces new contributions to the quantum corrections to the Higgs mass. The 
supersymmetry ensures that the contribution of new supersymmetric particles exactly cancel the 
quadratic divergence of their Standard Model partners, removing the sensitivity of the Higgs 
mass to the cut-off scale in a natural way. 

The Minimal Supersymmetric Standard Model (MSSM) contains a scale ^i, the Higgs-higgsino 
mass parameter in the superpotential, which is phenomenologically constrained to lie not far from 
the electroweak scale 3Ij in the Next-to-Minimal Supersymmetric Standard Model (NMSSM) this 
mass parameter may be linked dynamically to the electroweak scale in a natural way [5]-|lUj. 
The superpotential of the MSSM must be analytic in the fields, preventing the use of only one 
Higgs doublet for the generation of both up-type and down-type quark masses. Thus the model 
requires two Higgs doublets, which is also necessary to maintain an anomaly free theory. One of 
the doublets (Hu) provides a mass for up-type quarks while the other (-ff^) provides a mass for 
down-type quarks and charged leptons. The term ^HuHd in the superpotential of the MSSM 
mixes the two Higgs doublets. Since the parameter /i, present before the symmetry is broken, 
has the dimension of mass, one would naturally expect it to be either zero or the Planck scale 
(Mpi). However, if = then the form of the renormalization group equations implies that 
the mixing between Higgs doublets is not generated at any scale; the minimum of the Higgs 
potential occurs for {Hd) = 0, causing the down-type quarks and charged leptons to remain 
massless after symmetry breaking. In the opposite case, for /_i ~ Mpi, the Higgs scalars acquire 
a huge contribution to their squared masses and the fine tuning problem is reintroduced. 
Indeed, one finds that /i is required to be of the order of the electroweak scale in order to provide 
the correct pattern of electroweak symmetry breaking. 

The most elegant solution^ of this //-problem is to introduce a new singlet Higgs field 5, and 
replace the //-term ^H^Hd by an interaction \S{HuH(i). When the extra scalar field S acquires a 
^For other solutions to the /i-problem see Ref.|12|. 
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non-zero vacuum expectation value (VEV) an effective /i-term of the required size is automatically 
generated, and the effective // parameter may then naturally be expected to be of the electroweak 
scale, n = X{S). In this way fi can be linked dynamically to the scale of electroweak symmetry 
breaking (although the mechanism generating the common scale of order 100 GeV to 1 TeV 
is left unexplained). This situation naturally appears in the framework of superstring-inspired 

models ^Hj. At the string scale, the Eq symmetry can be broken, leading to the usual low 
energy gauge groups with additional U{1) factors (for instance SO{W) xU{l) x ...). These 
extra U{1) symmetries can be broken via the Higgs mechanism in which the scalar components 
of chiral supermultiplets acquire non-zero VEVs in such a way that all additional gauge bosons 
and exotic particles acquire huge masses and decouple, except for one singlet field S and a pair 
of Higgs doublets and Hd- 

However, this model still possesses an extra global U{1) symmetry, a Peccei-Quinn (PQ) 
symmetry |14| . which is explicitly broken in the MSSM by the /i-term itself. The breaking of this 
extra U{1) at the electroweak scale leads to the appearance of a massless^ CP-odd scalar in the 
Higgs boson spectrum, the PQ axion ^Sj. Unfortunately, this axion leads to astrophysical and 
cosmological constraints which rule out most of the allowed parameter space, leaving only a small 
window with 10"''' < A < 10"^° Due to the very small value of A, a very large value of (S) 
would be required to generate /i in the required energy range, and this model is unsatisfactory 
as a solution to the //-problem. However, this massless axion may be avoided by introducing a 
term cubic in the new singlet superfield in the superpotential. This term explicitly breaks the 
additional U{1) global symmetry, providing a mass to the CP-odd scalar. This model is known 
as the Next~to-Minimal Supersymmetric Standard Model (NMSSM), and is described by the 
superpotential 

W = u'' h^QHu - hdQHd - h,LHd + XSiHuHd) + ^kS"^- (1) 

Despite the removal of the PQ symmetry, the NMSSM superpotential is still invariant un- 
der a discrete Z3 symmetry: <I> — > e'^'^^^^^, where $ denotes the observable superfields. Left 
untamed, this Z3 symmetry would lead to the formation of domain walls in the early universe 
between regions which were causally disconnected during the period of electroweak symmetry 
breaking jl7j . Such domain wall structures of the vacuum create unacceptably large anisotropies 
in the cosmic microwave background jlSj. In an attempt to break the Z3 symmetry, operators 
suppressed by powers of the Planck scale can be introduced. However, it has been shown that 
these operators, in general, give rise to quadratically divergent tadpole contributions, which once 
again lead to a destabilization of the mass hierarchy \19^ . This problem can be circumvented by 
^The physical axion would acquire a small mass by mixing with the pion. 
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introducing new discrete symmetries to forbid or loop suppresses the dangerous tadpole contri- 
butions. In this case the breaking of the Z3 symmetry should be small enough to not upset the 
mass hierarchy but large enough to prevent the problematic domain walls |2U1 121j . Indeed, in 
these scenarios one would expect the surviving tadpole terms to be sufficiently suppressed to not 
effect the low-energy phenomenology described in this study. 

We analyze the mass spectrum of the NMSSM paying attention in this report not only 
to the lowest mass state but also to the heavier Higgs bosons. Characteristic mass patterns 
emerge for natural choices of the parameters. An important parameter is the extent to which 
the PQ symmetry is broken, measured by the size of the dimensionless coupling k. Moreover, 
the couplings of the Higgs bosons to the Z boson are studied for several representative scenarios 
appropriate to the production of the NMSSM Higgs bosons at the next generation of proton and 
electron-positron colliders. 

A central point of the paper is the analytical analyses of the Higgs mass spectrum and cou- 
plings in the approximation where the supersymmetry scale is large and the Higgs parameter 
tan /3 is moderate to large. The concise formulae which can be derived in this parameter range 
describe the system to surprisingly high accuracy while providing a valuable analytical under- 
standing of the NMSSM Higgs sector, which is significantly more complicated than the minimal 
extension. 

The paper is organized as follows. In Sec. (2) we describe the Higgs sector of supersymmet- 
ric models with an additional singlet superfield. The vacuum is examined and used to place 
constraints on the parameter space. The Higgs mass matrices are discussed, and an approxi- 
mate diagonalization of the CP-even mass matrix is used to provide analytic expressions for the 
masses and mixings, from which the couplings to the Z boson can be derived. These expressions 
are used to give a qualitative understanding of how the masses and couplings depend upon the 
parameters. In Sec. (3) the Higgs boson spectrum of the NMSSM is considered for three represen- 
tative parameter choices in detail, characterized by the breaking of the PQ symmetry. Natural 
values for the parameters are motivated and it is shown that the renormalization group flow from 
the GUT scale down to the electroweak scale prefers scenarios where the PQ symmetry is only 
slightly broken. In this case, two CP-even Higgs bosons and the lightest CP-odd Higgs boson 
have masses of order 100 GeV, while one CP-even, one CP-odd and the charged Higgs bosons are 
heavy. The couplings of the lightest Higgs bosons to the Z are reduced preventing their detection 
at LEP if kinematically possible otherwise. However, their observation is possible at a future 
linear collider operating in the TeV energy range with high luminosity, so that the NMSSM may 
be distinguished from the MSSM at such a facility. If the PQ symmetry is strongly broken for 
large values of k, the extra fields may become more massive and eventually decouple, making 
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such a comparison more difficult. The results are summarized in the Sec. (4). In the appendix, 
approximate solutions for the masses and couplings of CP-even Higgs bosons are deduced for 
physically appealing parameter domains. 

2 The NMSSM Higgs Sector 
2.1 The NMSSM Higgs Potential 

The NMSSM Higgs sector is distinguished from that for the MSSM by the addition of an extra 
complex scalar field, S. The Higgs fields of the model then consist of the usual two Higgs doublets 
together with this extra Higgs singlet. 

In the superpotential, already presented in Eqn.Q, the extra singlet is allowed to couple only 
to the Higgs doublets of the model, and consequently the couplings of the new fields to gauge 
bosons will only be manifest via their mixing with the other Higgs fields^. The superpotential 
leads to the tree-level Higgs potential 0: 

V = VF + VD + V,oiu (3) 

with 

Vf = \XSW\Hu\^ + \Hd\^) + \\HuHd + KS^\^ (4) 

Vd = Ifim'-M' + la'lHlH.W (5) 

Koft = mj, jHu\^ + mj,jHdf + ml\S\^ + [XAxSHuHd+^KA,S^ + li.c.], (6) 

where g = \/ g'^ + g''^ with g and g' being the gauge couplings of SU{2)l and U{V) interactions 
respectively, and adopting the notation HuH^ = eai3{Hu)'^{Hci)^ = H^H^ — H^H^. The first 
two terms, Vf and Vd, are the F and D terms derived from the superpotential in the usual way. 
We also include additional soft supersymmetry breaking terms, V^oft) by adding all the forms in 
the superpotential with arbitrary (dimensionful) coefficients^, and soft mass terms for the fields. 
This leads to the inclusion of the extra parameters: A\, Af^, mn^, "^J^d ^'^d rris- 

The parameter combinations XAx and may always be taken to be real and positive since 
their complex phases may be absorbed into global redefinitions of HuHd and S. However, A and 

■^Higgs self interactions and couplings to the sleptons, squarks and higgsinos will also gain extra contributions 

in the Lagrangian directly from the new terms in the Superpotential. 
^These definitions differ from those in Ref. 9 w.r.t. the sign of 
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K (and therefore Ax and rnay be complex. Significantly, if the combination Xn* is complex 
then CP will be violated at tree-level in the Higgs sector. This is in contrast to the MSSM where 
tree-level CP conservation is guaranteed by the structure of the Higgs potential. In this paper 
we will adopt the usual convention of assuming that all four parameters are real. For studies of 
the CP-violating NMSSM see Ref.|12]. 

While the essential new elements of the Higgs sector can be elaborated at the tree-level, one- 
loop contributions to the Higgs masses and couplings from top and stop loops are included in 
the quantitative analyses by introducing the term T^oop into the potential, with 



(Moop) 



32^2 





4 / 3 



2m, log -2 



Q2 2 



(7) 



The dependence on the Higgs fields is implicitly contained in the top, mt, and stop masses, ^, 
which are dependent, in total or partially, on the Higgs VEVs. For the stop masses, we use the 
tree-level relation, 
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^ - - ' m| + + 2mt ± J (m^ - m^)2 + ^.mjX^ ) , (8) 



*2/l ~ 2 V *^ V ^ Q ^' 

where we set mg = mu = MsuSY = 1 TeV, the running top quark mass mt = 165 GeV and, as 
an example, assume maximal mixing, Xt = ^/^MsusY- The two physical stop masses are given 
by mi_^ = 789 GeV and = 1196 GeV for these parameters. The renormalization scale, Q, is 
taken to be the running top quark mass. For the present analysis one-loop expressions for the 
masses and couplings suffice. The explicit expressions for the one-loop corrections to the mass 
matrices were adopted from Ref. j24j . 

The "MSSM limit" can be approached smoothly by letting A and k ^ on a linear trajectory 
(i.e. k/A constant), while keeping the /i parameter ^ = \vs/V2 and the parameters A\ (which 
plays the role of B in the MSSM) and fixed. In this limit the Higgs singlet field decouples 
from the system completely, and we regain the doublet Higgs sector of the MSSM in its canonical 
structure. 



The structure of the vacuum can be simplified by making a SU{2)l x U{1)y gauge transfor- 
mation to choose {H~) = 0. In order to ensure that (Hj) = is a suitable choice of vacuum, 
i.e. that the vacuum is uncharged, one also requires that it be a stable minimum^ with respect 
to small perturbations in , corresponding to a positive charged Higgs mass-squared: 

^The requirement ( ) = is trivially satisfied for (H^) — 0. 
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where Vu and Vd are the VEVs of -ff^ and respectively, multiphed by a factor of \/2. The 
resulting vacuum has been studied in Ref-IH], where it was shown that there are three degenerate 
vacua which differ in the complex phases of {S) and {H^). Here, we will restrict ourselves to a 
discussion of the vacuum where all the remaining VEVs are real, as described by, 

™ = 7l(o)' <^"> = 75C!)' <"> = 71°" 

with Vs, Vu, and Vd real and positive. 

For this vacuum to be a local minimum, we obtain three relations, linking the three soft mass 
parameters to the three VEVs of the Higgs fields: 

< = lf{vl-vj)-lx\l + kV2A, + .Vs)Xvs'-^-lx'vl (11) 
" 8 2 2 Vd 2 

= lgHvl-vl)-l>^^vl + liV2A^ + ^Vs)Xvs^-^X\l (12) 

rrig = -K Vg - -X V + nXvuVd + —^XAx -=KAf,Vs; (13) 

2 \'2 Vs \/2 



as usual, we have written v = yv^ + v^. The local stability of this vacuum is ensured by allow- 
ing only positive squared masses of the physical fields (i.e. the mass eigenstates) , which leads to 
useful constraints on the parameters of the potential. 

Proving this to be a global minimum of the vacuum is beyond the scope of this study. However, 
we consider all possible neutral CP-even vacuum states and ensure that the physical vacuum is 
of lower energy: 

(i) First of all, it is clear that (H^) = {H^) = {S) = is also a stationary point. It must 
be ensured that this symmetric minimum is not preferred to the physical vacuum with non-zero 
fields. This is possible in two ways. For the symmetric vacuum to be a stable minimum requires 
all of fn\^, ^'h^ m| to be positive^. Therefore, choosing the parameters such that one of the 
soft masses is negative, is sufficient. However, while sufficient, this constraint is not necessary. 
Alternatively, one may check that {V) < for the physical vacuum, making it a lower energy 
state than the symmetric vacuum for which {V) = 0. This constraint places an upper bound on 
the parameter A\. 

(ii) The only other neutral CP-even vacua which could be problematic are those where two of 

(i/^), (H^) and (S) are zero and the other non-zero. The requirement that the resulting vacuum 

^This constraint is different from the corresponding condition in the MSSM since (S) — 0, removing the effective 
/x-term. 
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energy {V) is greater for these points than for the physical vacuum can be fulfihed in the major 
part of the appropriate parameter space. Vacua where only one of {H^), {H^), {S) is zero lead 
to an over-constrained system, so that these vacua can only be realized for very specific choices 
of the parameters, and may be safely ignored. 



The NMSSM Higgs potential is automatically bounded from below for non-zero k. The two 
terms in Vp contain contributions which are quartic in the usual neutral Higgs fields, Hu, Hd, 
and in the new scalar, S, and will ensure that the potential is bounded from below. 

2.2 The Mass Matrices 

From the potential, the Higgs mass matrices and subsequently the mass eigenstates can be 
derived. After shifting the Higgs fields to the minimum of the potential (given by Ean. (|lU| )). 
they are rotated by an angle (3 in order to isolate the zero mass Goldstone states, G, which are 
absorbed by the W and Z bosons to provide their masses. 
For the charged fields these redefinitions can be written: 

HV = sin 3 — cos /3, , , 

(14) 

H+ = H+cos(3 + G+smp, 
where G~ = G'^* and = H^*. For the imaginary and real field components we have 

QmH^ = (Pisin/?-G°cos/3)/V2, 

9mi/0 = (Picos/3 + GOsin/3)/^/2, (15) 
9m S = P2/V2, 

and 

3fiei?° = (— 5isin 

^eHl = ( Sicosp + S2s\np + Vu)/V2, (16) 
^eS = {S^ + Vs)/V2, 
respectively. Choosing tan/3 = Vu/vd the zero-mass Goldstone modes decouple, and the resulting 
potential has terms for the non-zero mass modes given by 

( sA 

S2 ■ (17) 

The charged fields are already physical mass eigenstates with tree-level masses given by 
(!) : M|± =Ml + M^- \{\v)\ (18) 



^ma.. = Ml±H+H- + ^ (Pi P2) M2 f + 1 (5i 52 ^3) 
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where 

^'^ = ii^(^-*^ + <'^' 

In contrast, the pseudoscalar and scalar fields, Pi {i = 1, 2) and Si {i = 1, 2, 3) respectively, 
are not yet mass eigenstates. Their mass matrices, and M^, must be further rotated to a 
diagonal basis corresponding to the physical mass eigenstates. Using the minimization conditions, 
Eans. lfTTHT^ . the tree-level CP-odd matrix has entries, 

(//): = Ml, (20) 

M^i2 = ]^(M\sm2(3 -^\Kvl)coi(3s, (21) 
M^ga = ^(M|sin2/? + 3Araf)cot2/?sSin2/?-3CT^A/\/2. (22) 

while the tree-level CP-even matrix is given by, 

(///): Ml^^ = Ml + (M|-i(At;)2)sin2 2/3, (23) 

Mli2 = -\{Ml-\[\vf)smA(5, (24) 

M^i3 = -^(Mlsin2/3 + ACT^)cot/?^cos2/3, (25) 

M^22 = Mlcos^2[3 + ]^{\vfs\n^2l3, (26) 

M^23 = ]^{2\^vl - M\s\ii^ 2(5 - XkvI sm2(3) cot 13s, (27) 

1 1 

M^33 = -M\sm^ 2(3 cot^ 13s + 2k'^vI + KVsA^/V2- -\Kv'^sm2l3, (28) 

Besides the the usual notation for tan/3, we have introduced the auxiliary ratio tan/J^ = Vs/v. 
The mass parameter M\ is defined to be the top-left entry of the CP-odd squared mass matrix, 
c.f. Ean. ()2U( ). which is positive in the physical parameter ranges analyzed later. It becomes the 
mass of the heavy pseudoscalar Higgs boson in the "MSSM limit". Parallel to the MSSM, Ma 
replaces the soft parameter A\, and at tree-level is given by Ean. lfT^ . 



Higgs Mass Spectrum 

The charged Higgs mass has been noted in Eqn. (|18)) . The condition for stability of the vacuum 
in the Hj direction, Eqn.Q, is simply equivalent to the positivity of M'^±. 

The above pseudoscalar and scalar mass matrices, Eqns. (|20H28() . do not lend themselves 
easily to obtaining analytic expressions for the physical Higgs masses. However, reasonably simple 
expressions can be found by performing a systematic expansion for large tan /3 and large M^ [on 
the generic electroweak scale]. This is outlined in detail in the appendix. This approximation 
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works extremely well, even for moderate values of tan /?, and it may therefore be used to shed 
light on the behaviour of the Higgs masses as the other parameters are varied. 

This approximation may be used to simplify the expressions for the tree-level CP-odd masses , 
giving, 

Ml « Ml (l + ^cot2/?,sin2 2/3), (29) 
Ml ^ -^KVsA^, (30) 

while the tree-level physical masses of the CP-even Higgs bosons are, in this approximation, 

1 
4 



M^3 « Ml(l + -cot2/3,sin2 2/3), (31) 



If 1 

Ml^^^ ^ -lMl + -KVs{^KVs + V2A^) 



±1 



+ cot2 /3, [2A2t>2 - M\ sin2 2/3] ^ J. . (32) 



The physical mass eigenstates Ai and Hi are labelled in ascending order of mass. 

The heavy CP-odd Higgs boson, A2, is approximately degenerate with the heaviest CP-even 
Higgs boson, i/3, and the charged Higgs bosons, cf. Ean. H18|) . One of the lighter CP-even tree- 
level masses will be of the order of Mz while the scale of the other is set by ~ kVs for k and 
A sufficiently below unity. Finally the lightest pseudoscalar mass grows as k, Vs and A^ are 
increased, with a negative value for A^, being preferred. 



Range of Parameters 

This above solution also allows limits to be placed on the Higgs potential parameters. First 
of all, the lightest two scalar Higgs boson masses respect a sum rule: 

Ml + Ml ^Ml + ]^KVs{^KVs + V2A^). (33) 

The right-hand side, and thus the sum of the two lightest scalars, is independent of the coupling 
A and M^. As a result of this sum rule, the second lightest scalar Higgs boson mass is maximized 
as the lightest approaches zero. 

Furthermore, the smallest eigenvalue of a matrix is smaller than its smallest diagonal entry 
and similarly the largest eigenvalue is larger than its largest diagonal entry. Together with 
the above sum rule, this leads to the mass constraints at tree-level (i.e. modified by radiative 
corrections in parallel to the MSSM): 

Ml <mm{Ml]^KVs{^KVs + V2A^)}, (34) 
max{/Vf|,iCT,(4CT, + ^/2A«)} < Ml < Ml + ]^kVs{^kVs + V2A^). (35) 
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A further constraint is found by exploiting the condition that the mass-squared of the hghtest 
Higgs boson, given in Ean. (|32| ). is greater than zero. This in turn gives a restriction on the allowed 
values of Ma- At tree-level, again to the accuracy of the approximate solution, this is, 

Ml S ^ ± (4-. + V/2A.), (36) 

sm Ifj sm 2p \ 2 

which generalises the corresponding constraint for zero n in Ref. In the theoretically pre- 
ferred scenarios where k is small, the second term above is small compared to the first term, and 
we find that the value of M^ is constrained to lie in a narrow bracket, not too far from the value 
■\/2Xvs/ sin 2/3 = 2fi/ sin 2/? which is approximately ~ /utan/3 for medium to large values of 
tan/?. As can be seen from Ean. (|32j l. this is also the approximate value of Ma for which Mhi 
is maximal (and Mh2 minimal). The constraint on Ma becomes stronger as k becomes smaller 
but Vs is kept fixed. The mass parameter Ma is not constrained in the proper "MSSM limit" 
defined earlier with Xvs and kVs fixed for A and n ^ 0. 

Finally, one may also gain some insight into the allowed range of the soft SUSY breaking 
parameter A^- Requiring the lightest scalar and pseudoscalar mass-squareds be positive leads to 
upper and lower constraints on A,^. At tree-level, these are given by, 

- 2V2^,Vs < A < 0. (37) 

The lower bound is derived from the requirement that the maximum value of the lightest scalar 
mass-squared be positive, which is realized for the central value of Ma given in Ean.H36|l. For 
values of Ma deviating from this, the lower bound will become stronger. 

One should bear in mind that most of these constraints are based upon the approximate 
solution as described in the appendix, and will become unreliable in regions of parameter space 
where the approximate solution breaks down. Nevertheless they provide the proper analytical 
understanding of the exact numerical results presented in Sec. (3). 

2.3 Couplings with the Z boson 

The couplings of the Higgs bosons to the Z boson are given by the Lagrangian, 

CzAH = ^MzZ^Z''S2 + ^Z^[Si{d^Pi)-Pi{d''Si)] 

+ (eA^ + cos 201^) i [H-{d^'H+) - H+{d^H-)] , (38) 

where Z^ and are the Z boson and photon fields respectively, 9w is the weak mixing angle 
and the Higgs fields have been defined in Eans. ()14H16() . The simple form of this Lagrangian 
arises from the rotation made on the CP-even Higgs states. This leads to 5*2 being the only field 
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with a coupling of the form ZZSi, while the scalar and pseudoscalar fields are coupled jointly to 
Z only in the form ZS\Pi. The couplings of the extra scalar and pseudoscalar Higgs bosons to 
the Z are due only to their mixing with the Higgs doublet degrees of freedom. We define this 
mixing of the CP-odd Higgs bosons via a rotation by an angle 9a which transforms the fields to 
their physical mass eigenstates: 




cos 6 a sin Oa 
— sin 6 a cos 6 a 




(39) 



At tree-level, the mixing angle, 6a, is given by. 



tan 6'yi 



12 



1 ^ M'\ sill 2, -i - 3Ak;'2 

— cot Do — 7^ r; . 

2 Ml- Ml 



(40) 



Similarly, the more complicated scalar mixing is defined via an orthogonal 3x3 rotation 
matrix O such that, 

f ( Si\ 

= O 



H2 



S2 



(41) 



Our notation is such that the fields Ai, i = 1,2, and H^, i = 1,2,3 are the mass eigenstates of the 
pseudoscalar and scalar Higgs sectors respectively, with the numerical suffix denoting their mass 
hierarchy in ascending order (e.g. Ai (Hi) is the lightest pseudoscalar (scalar) Higgs boson). 
[The ordering of the states Si has been introduced such that S3 is the decoupling singlet.] 

It is more useful to discuss the couplings of these Higgs bosons to the Z in terms of normalized 
couplings. To this end, we define the normalized CP-even couplings by dividing out the associated 
SM coupling and for the ZAiHj couplings we similarly divide by g/2 = iyp~+~3'^/2, so that. 



QzzHi = 



„NMSSM 

9zzH 



QzAiHi = 



„NMSSM 
y ZAiHj 

5/2 



(42) 



The couplings may then be written directly in terms of the mixing matrices, in an obvious 
notation, according to, 

Gzzh, = Oh,s„ (43) 

and 

QzAiHi = -OHiSi sine A, QzA2Hi = OuiSi COS^A- (44) 
For these normalized couplings the orthogonality of the mixing matrices lead to the constraints, 

^QzZHi = 1> ^QzAiHj = 1> ^QzZHiQzAjHi = 0. (45) 
i ij i 
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Finally, the couplings of a pseudoscalar and a scalar to a Z boson are not all independent. The 
ratios of the two pseudoscalar couplings, Ai and A2, are independent of the scalar Hf. 

QzAiHi _ QZA1H2 _ GzAiHs 



GZA2HI GZA2H2 GZA2H3 



tan 9 A. (46) 



In the approximate solution, where 1/Ma and l/tan/3 are regarded as small quantities, the 
orthogonal matrix O is given by Eqn.ljTHJ in the appendix. Inserting this into Ean. (|43() . gives 
approximate expressions for the scalar- Z-Z couplings. 



GzzHi 



sm t7H, 



QzzH2 ~ COS Oh-, 



where tanOH is given by (see Ean. ((7^ ). 



2M\ 



\{Xvf-Ml 



sin 4/3, 



1 ^ 2yvt - Mi sin^ 2/3 - nXvi sin 2/3 
tanOu - cot Ps — K 7^ 1 — TT-^ n 7^- 

2 M| cos2 2/3 + ^X^v^ sin2 2/? - M^^ 

Similarly, for the Z-pseudoscalar-scalar couplings we have, 

I 

GzAiHi ~ [M\ sin 2(3 + kXvI] cot (is cos 2/3 sin 61a, 



GzAiH2 
GzAiH-i 



2M\ 



1 + 



-iXvf-Ml 



sin 4/? sin 0^, 



4M| 



[Ml sin 2/3 + kAz;^] cot (5s cos 2/3 



sin 



and, following from the relations Ean. (|46( ) 

1 



GZA2H1 
GZA2H2 
GZA2H3 



4Ml 



[M\ sin 2(5 + kXvI] cot (5s cos 2/3 cos 9a, 



1 



2M| 



(At;; 



|2 - m| 



sin 4/3 cos 9 a-, 



1 + [Ma sin 2/3 + kA^;^] cot (5s cos 2/3 



cos 9 a, 



where tan^^ is given by Ean.H40|). 



3 NMSSM Higgs boson scenarios 



(47) 
(48) 

(49) 

(50) 
(51) 
(52) 

(53) 
(54) 
(55) 



3.1 The parameters 

At tree-level the NMSSM Higgs sector described above has six free parameters (apart from the 
overall scale, set later by the Z boson mass). The Higgs potential, Eqns.(I3Hni)i contained seven 
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parameters: A and k from the superpotential and A^, mn^, itt-Hu ^''^^ from the soft 
supersymmetry breaking terms. The field values at the minimum of the potential, Vu, Vd and Vg, 
are fixed by these parameters according to Eans. (|llH13() . The structure of electroweak symmetry 
breaking allows us to remove one (combination) of these VEVs in favour of the known electroweak 
scale V = y''^2 + = 246 GeV (or equivalently, the Z boson mass, Mz), defining the overall 
mass scale. Finally, after introducing tan [3 = Vu/vd and writing Ax in terms of the heavy Higgs 
mass scale Ma, we are therefore left with the six free parameters, A, k, Vs, tan/3. Ma, and 
{vs may be expressed in some places in terms of the effective /i parameter, Vg = V^fi/X). The 
spectrum of the NMSSM is expected to show a strong dependence on these six free parameters, 
which we shall now analyze in detail. 

Including higher orders introduces extra parameters such as the top and stop (s)quark masses. 
We choose to fix these extra parameters at reasonable values and not vary them. 



A and k 



Requiring a weak coupling of the fundamental fields, i.e. field-theoretic perturbativity, in 
the entire range between the electroweak and GUT scales restricts the range of values for the 
couplings A and k at the electroweak scale. The renormalization group equations for A, k and 
the top Yukawa coupling ht form a closed set together with those of the gauge couplings. They 
are given by El 



16vr -TT = Mi, 



167r 



dt 
2 dhl 
dt 
d>?_ 
dt 

dK^ 



A^ + 6/.?-^5i-3,|-H,2 



A^ 



4A2 + 2k^ + 3/i? 



35i 



3 2 
5'^ 



dt 



(56) 
(57) 

(58) 

(59) 



[A^ + K^] , 

where bi = 33/5, 62 = 1, 51 = \/5/3g', 52 = 5 and t = log{Q^ / M^^.^) . 

Large values of A and/or k at the GUT scale are greatly reduced when run down to the 
electroweak scale. This behaviour is caused by the dependence on A'^ and on the right-hand 
side of Eans. (|58H59|) respectively, indicating that large values of these parameters will evolve 
strongly, while small values evolve only slightly. This can be seen in Fig.(^, which shows the 
dependence of A and k. on renormalization scale. Values of A and k in the perturbative regime 
at the GUT scale, i.e. A, k < 27r, are uniformly reduced to small values at the electroweak scale, 
which may be combined to give the approximate bound, c.f. Fig.(|2/left), 

X^ + K^< 0.5. (60) 
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Figure 1: The dependence of A (left) and k (right) on renormalization scale, Q, for the top 
Yukawa coupling /i((Mgut) = 0.8. The different curves represent different values of X and k at 
the GUT scale; left: A = 1, 2 and 3, with k = 1; right: k = 1, 2 and 3, with A = 1. 

Also notice the factor of 6 in Ean. (|59( ). compared to the smaller factors in Ean.()58(). This tends 
to make k run more steeply than A, which is further exacerbated by the lack of gauge couplings 
on the right-hand side of Eqn.(|59j). since the k term in the superpotential contains only the 
superfield S which has no gauge couplings. 

In order to illustrate these behaviours more clearly, we considered 2 x 10^ different parameter 
scenarios with (positive) GUT scale values of A, k and ht chosen randomly (and independently) 
between and 27r, corresponding to the perturbative regime defined by k/2'k < 1, A/27r < 1. 
Using the renormalization group equations, Eans. (|56H59|) . to run A and k down to the electroweak 
scale, the final distribution in A and k can be seen in Fig.((21/right), where each point is a different 
parameter choice at the GUT scale. Only scenarios where the running top quark mass mt{Q) 
falls in the bracket 165 it 5 GeV at the electroweak scale are retained (approx. 6%). This plot 
demonstrates the limit on A^ + k^ of Ean. (|6U| ) since no parameter point strays above A^ + k^ ~ 0.6. 
Furthermore, one can easily see that most of the parameter points lie at low k. Nevertheless, the 
entire area within A'^ + < 0.6 is valuable; the accumulation at the boundary results from the 
mapping of k and A from the large number of moderately high values at the GUT scale, to small 
values at the electroweak scale. If n is small at the GUT scale then its renormalization group 
running will be weak and it will remain small at the electroweak scale. 

The size of k governs how strongly the U(l) PQ symmetry of the model is broken. From the 
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Figure 2: Left: The dependence of + on renormalization scale, Q, for the top Yukawa 
coupling /i((Mgut) = 0.8. The different curves represent different values of X and k at the GUT 
scale: X = k = 1, 2, and 3. Right: The electroweak scale values of k and X for 2 x 10^ different 
scenarios with random GUT scale values o/O < (A, n, ht) < 27r. Only scenarios where the running 
top quark mass mt{Q) falls in the bracket 165 it 5 GeV at the electroweak scale are retained. Each 
point represents a different GUT scale parameter choice. 

above investigation, it is natural to expect that small values of k are preferred and this symmetry 
is only 'slightly' broken. However, this is not a strict bound; scenarios with k 0.6 can still 
occur. 

However, if one wants the NMSSM to be a solution to the /i-problem, A must not become too 
small. We have introduced Vs as the VEV of a new Higgs field in order to link to the electroweak 
(or supersymmetry) scale. If A becomes too small then the phenomenological constraints on ^ 
require Vs to be large and the link with the other Higgs VEVs is lost. Even if we allow Vg as 
large as 3 TeV (which is a natural upper bound for a weak scale) we would be forced to maintain 
A > 0.04. Consequently the allowed values of A are rather constrained. 

Even more stringent correlations among the parameters are imposed by grand unified scenar- 
ios of the NMSSM including universal boundary conditions for the soft masses and the trilinear 
couplings. Radiative electroweak symmetry breaking requires the couplings, at the unification 
scale, to be less than about A^, and both less than about 0.1 [7]. To prevent QED and color 
breaking vacua, the ratio of the universal trilinear coupling and the scalar mass parameter is 
restricted to values close to 3 0. If universal boundary conditions are not imposed, these con- 
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straints will not be effective, of course. 
Vs and tan P 

The Higgs-higgsino mass parameter /i is generally assmned to be in the range ~ 10^-10^ GeV. 
However, the 'natural' scale for the heavy Higgs masses Ma is seen by Eqn. (|.S6) ) to be of the order 
of n tan p. Indeed, this is approximately the value of required for the lightest scalar to reach 
its maximum mass. In order to prevent the mass splitting between the light and heavy Higgs 
bosons from becoming too large^, the value of tan/3 should be kept moderate, < 10. The value 
of Vs could naturally be expected to vary from 3f to 15 v. 

A similar analysis of the renormalization group equations for the other parameters shows 
that a low value of tan /? is favoured 26 . Note that the bounds on tan /? derived from the LEP 
experiments cannot be applied here since these experiments assume MSSM couplings. Referring 
to the literature for the large tan P case |27| . here we will adopt a value of tan /3 = 3. 

and Ma 

The value of A,^ is tightly constrained by the requirement for vacuum stability, as described 
in Sec. (2.2). Note that it occurs in only one term of each of the CP-even and CP-odd mass 
matrices, Eqns. lj^l^ . where the two contributions have opposite sign. U A^ becomes too large 
and positive, it will pull the mass-squared of the lightest pseudoscalar below zero, destabilizing 
the vacuum. On the other hand, if it becomes too large and negative it will destabilize the vacuum 
by pulling the mass-squared of the lightest scalar negative. This has been quantified in Ean.()37|l 
and can be seen in Fig.© where the Higgs masses are plotted as a function of for reasonable 
values of the other parameters, including radiative corrections. The heavy Higgs bosons and 
one of the lighter scalars are very insensitive to the choice of Af^, but the predominantly singlet 
scalar and pseudoscalar Higgs bosons exhibit a strong dependence. This is naturally explained 
by noting that the soft SUSY breaking term containing is a cubic coupling proportional to , 
so the affect of varying is only communicated through the singlet contribution of the Higgs 
fields. Since the mixing is small, this only affects the two fields which are still singlet dominated. 

In contrast to the MSSM, Ma is also constrained by vacuum stability in the NMSSM, and, 

as discussed above and visualized in later plots, its "natural" value is approximately /itan/3. 

Nevertheless, we will allow Ma to vary over the small bracket of its allowed range, including the 

"natural" value, and plot the Higgs masses as a function of Ma- 

^Indeed, if a fairly large splitting were to be observed in experiment, then the NMSSM with large fi and/or 
tan/3 would provide a very natural explanation. 
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Figure 3: The one-loop Higgs boson masses as a function of for A = 0.3, n = 0.1, Vg = 2>v, 
tan/? = 3 and Ma = /itan/3 f» 470 GeV. The arrows denote the physically allowed region. 

3.2 The NMSSM with the Peccei-Quinn U(l) Symmetry 

For illustrative purposes, first consider the simplified case where the U(l) PQ symmetry is left 
unbroken in the superpotential, achieved by setting k = 0. The U(l) symmetry is then sponta- 
neously broken by the vacuum, giving rise to a massless Goldstone boson, the PQ axion, which is 
manifest as the extra pseudoscalar Higgs field. Since there is no SUSY breaking term proportional 
to Afi, there are only 4 parameters in the model at tree-level: A, Vg, tan/3 and Ma- 
While the charged Higgs boson masses are, as usual, given by 

(/): Ml^=Ml + M^-]^{\v)\ (61) 

the tree-level pseudoscalar Higgs boson masses read, 

(//): Ml^ = Ml (l + ^cot2/?,sin2 2/3), (62) 

Ml, = 0, (63) 

where the expressions are exact, with no need of approximation. 

The approximate solution, where 1/M^ and 1/ tan [3 are regarded as small parameters of 0(e), 
may again be used to shed light onto the behaviour of the CP-even masses as the parameters 
are varied. However, due to accidental cancellations, the approximate solution taken to the same 
order as Ean. H32( ). O(e^), always yields a lightest scalar mass-squared that is less than zero. To 
circumvent this, one must improve the approximate solution to order O(e^). For the scalar Higgs 
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masses, we find 



(III) : ^ Ml (1 + i cot^ /3, cos^ 2/3 sin^ 2/3) + (M| - X^v^) sin^ 2/?, 



(64) 



M 



Ha 



M 



M| + Jm| + cot2 /3, [2A2?j2 



Ml sin2 2/?] ■ 



(Ml- AV)sin2 2/3,(65) 



M| - a/ M| + cot2 /?, [2A2t;2 - Ml sin2 2/?] ' 



+ -X^v"^ sin^ 2/?, 



(66) 



where in the last term of each, i.e. the terms of 0{e'^), we have made the replacement 
Ma = \/2\vs/ sin 2(3 for simplicity of the expressions. These relations improve on those of Ref. 
by describing the characteristics of the Higgs mass spectrum away from the mid point of Eqn. (|36| ). 

The one-loop Higgs boson masses are plotted as a function of Ma with A = 0.3, k = 0, = 
and tan/3 = 3 in Fig.Q, including radiative corrections. The five heavy Higgs bosons, H^, H2, 




Figure 4: The one-loop Higgs boson masses, plotted as a function of Ma for A = 0.3, k = 0, 
Vs = 3v and tan /? = 3. The arrows denote the physically allowed region. 

A2 and H^, show a pattern very similar to the MSSM. However the spectrum is augmented by an 
additional massless pseudoscalar Goldstone boson and a very light scalar Higgs boson, conforming 
with the tree-level upper bound, Ausin2/3/\/2, as evident from Eqn. (|66|) . The positivity of the 
mass-squared of this additional scalar Higgs boson places a very stringent bound on the value of 
Ma, predicting a small bracket in analogy to Eqn.(|36j) for the k non-zero case, and consequently 
strongly constraining the masses of the Higgs bosons H^, A2 and H^. 
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Taking A — > 0, besides the massless pseudoscalar PQ axion, also the mass of the additional 
scalar Higgs boson approaches zero. Both fields decouple from the doublet Higgs sector. At the 
same time, however, the degenerate masses of the heavy Higgs bosons are fixed to the approximate 
value Ma ~ /itan/5. For k = 0, the "MSSM limit" is approached on a trajectory A ^ that is 
constrained to a sheet in the MSSM parameter space that respects the PQ symmetry in the final 
MSSM potential. In fact, the HuH^ mixing term, that violates the PQ symmetry in the Higgs 
potential, is suppressed compared to the diagonal term, which respects the PQ symmetry, 

in the limit of large Ma_ and tan/3 analyzed analytically above. 

While the present model may serve as a nice illustration for the complex phenomena in the 
NMSSM, the presence of a massless pseudoscalar and light scalar rules out this parameter choice, 
due to their couplings to the Z boson. The direct production of the lightest scalar together with 
a Z boson, Z ZHi, the coupling of which is shown in Fig.((2^'left), would have been possible at 




Figure 5: The couplings QzzHi Oeft) and QzAiHj (i^ight), plotted as a function of Ma for A = 0.3, 
K = 0, Vs = 3v and tan/? = 3. The arrows denote the physically allowed region. 

LEP for much of the Ma range. However, since the coupling passes through zero, such production 
alone does not yet rule out the model. The Z boson decay Z HiAi is more damaging, since the 
coupling QzAiHi (see Fig.ljnj^right)), although small, is always large enough to allow detection. 

This may only be remedied by forcing A to become very small, so that the extra scalar 
and pseudoscalar fields decouple from the SM gauge bosons. Disregarding potential fine-tuning 
problems, these models are further constrained by astrophysical axion searches and cosmological 
bounds which are only satisfied for IQ-^o < A < 10"^ Since /i is constrained to be of the 
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order of the electroweak scale, this requires the extra singlet field to have a very large VEV, 
making this model unattractive as a solution to the ^u-problem. 

3.3 The NMSSM with a slightly broken Peccei— Quinn symmetry 

Turning on a non-zero value of k breaks the PQ U{1) symmetry, providing a non-zero mass for 
the lightest pseudoscalar Higgs boson, and raising the mass of the lightest scalar Higgs boson. 
We consider the PQ symmetry to be only 'slightly' broken as long as k ^ 1, for values of Vg of the 
order of the electroweak scale. This is favoured by the renormalization group flow as illustrated 
in our earlier discussion. The model contains two extra parameters as compared to the case with 
an unbroken PQ symmetry: k and its associated soft SUSY breaking parameter A^. 

As a representative example of this scenario we set A = 0.3 and k = 0.1. We will choose a 
value of in the centre of the allowed range, A^^ = —100 GeV, c.f. Fig.©, but the masses of 
the two singlet dominated fields can change significantly by altering this value. 

The masses of the Higgs bosons for A = 0.3, k = 0.1, Vg = 2>v, tan/3 = 3 and = —100 
GeV are shown as a function of Ma in Fig.®, including radiative corrections. The structure of 




Figure 6: The one-loop Higgs boson masses as a function of Ma for A = 0.3, k = 0.1, Vg = 3v, 
tan/3 = 3 and A^ = —100 GeV. The arrows denote the region allowed by LEP searches with 95% 
confidence. 

the heavy Higgs boson masses, i?3, A2 and H^, remains very similar to the MSSM. However, 
the light Higgs boson mass spectrum is distinctly different, with two light scalar Higgs bosons 
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and one light pseudoscalar Higgs boson. The extra singlet dominated pseudoscalar is no longer 
massless. Its mass has been raised from zero, by having the PQ symmetry explicitly broken by 
the term ^kS^ in the superpotential, to a value largely independent of Ma, c.f. Ec pUf) . Once 
again, the value of Ma is bounded by the requirement that the vacuum be stable (i.e. > 0), 
although the restriction is now looser. The H2 mass, of the order of Mz and largely independent 
of Ma, can rise above the canonical MSSM value due to the interaction with the new singlet 
fields, c.f. Ea (P^ . Generating a spectrum of light Higgs states, the NMSSM would be easily 
distinguished from the MSSM for this parameter range except in small regions where the ZZHi 
couplings are particularly close to zero. 

Also shown are the restrictions on Ma imposed by the LEP experiments [^Hl with 95% confi- 
dence. These limits are dependent on the masses of the light Higgs bosons and their couplings to 
the Z, as shown in Fig.((7J. The light scalar escapes detection via Higgs-strahlung only when it 
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Figure 7: The couplings GzzHi (left) o,nd GzAiHj (right), as a function of Ma for A = 0.3, 
K = 0.1, Vs = 'iv, tan/3 = 3 and = —100 GeV. The arrows denote the region allowed by LEP 
searches with 95% confidence. 

has a sufficiently reduced coupling to the Z boson. Since this coupling passes through zero in the 
allowed range, one is unable to rule out all values of Ma- Furthermore, in contrast to the k = 
scenario, the lightest pseudoscalar is now heavy enough to escape production via e'^e~ HiAi 
at LEP. If the couplings of the lightest Higgs boson are particularly small, this could be difficult 
to discover even at the LHC^. It is therefore important that a future linear collider, such as 
*The couplings of the new fields to quarks and leptons are manifest only from their mixing with the usual Higgs 
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JLC/NLC/TESLA [SOl-ISl] or CLIC |33] perform a search for light Higgs bosons with reduced 
couphngs^. 

Varying the parameters Vg and changes the quantitative, if not quahtative, features of 
the mass spectrum. In particular, large departures from the typical MSSM Higgs mass pattern 
occur in the light sector. As pointed out in Sec. (3.1), changing the value of varies the lightest 
pseudoscalar and scalar masses while leaving the heavier Higgs bosons unaffected (see Fig. 

If Vs (or A) is increased, all but the lightest Higgs boson become heavier. This can be seen 
in Fig.(I5}^left) where the Higgs boson masses have been plotted with Vg = 10 v, A^ = —250 GeV 
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Figure 8: The one-loop Higgs boson masses(left) and couplings QzzHi (fight) as a function of 
Ma for A = 0.3, K = 0.1, Vg = lOw, tan/3 = 3 and A^ = —250 GeV. The arrows denote the 
region allowed by LEP searches with 95% confidence. 

and the other parameters as before^'^. The restrictions on Ma enforced by vacuum stability 
(i.e. > 0) are relaxed, allowing a larger range of Ma values. One should note that 

Ma ^ y/2\vs/sm2f5 « ^tan/3 remains a point of particular interest as this choice allows the 
largest mass for Hi. Ma is now forced to be beyond the 1 TeV range. 

doublet fields, just as for their couplings to the gauge bosons. Small mixings will therefore also suppress the (loop) 

coupling of the predominantly singlet Higgs bosons to gluons. Searches for Higgs decays to light axion pairs at 

colliders have been discussed in Ref.|29). 

^For the phenomenology of Higgs-strahlung, e^e~ ZHi, and associated scalar/pseudoscalar production, 

e^e~ — > HiAj, for heavy and light Higgs bosons, see Ref.l34|. 
^^Ak has been chosen to lie in the middle of its allowed range. 
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Significantly, the increase of Vg causes the mixing to become stronger, increasing the doublet 
component in Hi at the expense of the doublet component in H2- Consequently, the character 
of Hi becomes more doublet-like, while the character of H2 becomes more singlet-like. This in- 
creased mixing is quantified by Ean. (|in)) . Consequently, the lightest scalar Higgs boson acquires 
a significant coupling to the Z boson, as shown in Fig.(jSl/right). 

A further example of models in this class is given by A = 0.05, k = 0.02. Such a small value 
of A requires a large value of Vg, here taken to be 15?;, in order to satisfy the phenomenological 
constraints on fi. The Higgs mass spectrum and couplings for this model remain qualitatively 
the same, as can be seen in Fig.®, although a larger range of Ma values remains viable. 




(GeV) (GeV) 

Figure 9: The one-loop Higgs boson masses (left) and couplings GzzHi (right) as a function of 
Ma for A = 0.05, k = 0.02, Vg = 15 v, tan/3 = 3 and = —100 GeV. The arrows denote the 
region allowed by LEP searches with 95% confidence. 

Finally, changing the sign of k, while keeping that of kA^ fixed, results in behaviours very 
similar to their positive k counterparts. Changing the sign of while keeping that of k fixed 
results in an unstable vacuum, as can be seen in Fig.©. 

3.4 The NMSSM with a strongly broken Peccei— Quinn symmetry 

When K becomes large the PQ symmetry is strongly broken. This scenario has been shown earlier 
to be disfavoured by the renormalization group flow but cannot be ruled out a priori on general 
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grounds. The parameters Vg, Ma etc. are taken to be of the order of the electroweak symmetry 
breaking scale. The extra pseudoscalar and scalar Higgs bosons gain moderately large masses 
from the PQ breaking term ^kS^ in the superpotential, with values M^^ ~ — Sku^^k/v^ and 
~ kVs{kVs + \/2j4k)/2, c.f. Eans. (|5n|) and ((32), largely independent of Ma- The eigenstate 
predominantly composed of the new scalar may in general no longer be the lightest scalar Higgs 
boson. That role is taken over by a Higgs boson which is predominantly composed of the doublet 
field 5*2, with a mass of order Mhi ~ Mz and a SM like coupling to the Z. The masses of these 
Higgs bosons are shown as a function of Ma for the parameter choice A = 0.3, k = 0.5, Vg = 3v, 
tan/3 = 3 and Af^ = —500 GeV, in Fig. (|10() . Again the predominantly doublet fields follow the 
pattern of the MSSM, with one CP-even Higgs boson. Hi, below the electroweak scale and one 
CP-even, one CP-odd and two charged Higgs bosons nearly degenerate at the scale Ma- However, 
now the predominantly singlet scalar and pseudoscalar fields are heavy enough that they are no 
longer the lightest states for much of the physically allowed range, being, independently of Ma, 
at the scale of a few hundred GeV. 

The couplings to the Z boson. Fig. 1)11(1 . show a complicated structure, due mainly to the 
way in which the Higgs states are labeled. As AIa is increased a previously light Higgs boson 
may gain sufficient mass to overtake one of the other Higgs bosons. The labeling will change 
(according to the mass hierarchy) causing an apparent rapid variation of the Z couplings at the 
crossing. It is this relabeling which gives the complicated structure of the ZAiHj couplings. 

Increasing the value of Vs will increase the mass of the new fields substantially (since the effect 
of the increase seen in the previous section is magnified by the large value of k). For example, 
with the same parameters as Fig. (|10|) except Vg = 10 v, the new scalar and pseudoscalar states 
have masses of the order of 1.6 TeV and 1.2 TeV, respectively, for all values of Ma (although 
these masses show a strong dependence on Af^). The new fields will also decouple, making the 
NMSSM very difficult to distinguish from the MSSM. 

Once again, negative values of k do not change the qualitative picture as long as one also 
changes the sign of A^- 

4 Summary and Conclusions 

In this study, we have investigated the Higgs sector of the Next-to-Minimal Supersymmetric 
Standard Model, suggested by many GUT and superstring models. Moreover, this model at- 
tempts to explain the //-problem of the MSSM by introducing a new singlet Higgs field, S, with 
a non-zero vacuum expectation value. When this new field is coupled to the Higgs doublets, its 
expectation value leads to the /t term of the MSSM, providing a value of /x which may naturally 
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Figure 10: The one-loop Higgs boson masses as a function of Ma for the parameter choice 
A = 0.3, K = 0.5, Vs = 3v, tan/3 = 3 and = —500 GeV. The arrows denote the region allowed 
by LEP searches with 95% confidence. 
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Figure 11: The couplings QzzHi (l^ft) and QzAiHj (i^ight), as a function of Ma for the parameter 
choice A = 0.3, k = 0.5, Vs = 3v, tan/? = 3 and = —500 GeV. The arrows denote the region 
allowed by LEP searches with 95% confidence. 



25 



be expected near the electroweak scale. 

We have given expressions for the Higgs boson mass matrices and presented, besides the 
numerical analyses, approximate analytical solutions for the charged, CP-even and CP-odd Higgs 
boson masses which provide a nice insight into the mass hierarchies. It has been found that useful 
sum rules and inequalities for the mass parameters can be established and that the requirement of 
vacuum stability can place useful bounds on the mass parameters and couplings of the Lagrangian. 

The renormalization group flow of the parameters A and k from the GUT scale down to the 
electroweak scale provide strong upper bounds on their values at the electroweak scale, where 
small K is favoured. The qualitative features of the Higgs boson masses are dependent on how 
strongly the PQ symmetry of the model is broken, quite accurately described by the approximate 
analytical solutions. 

If the PQ symmetry is left explicitly unbroken in the Lagrangian (spontaneously broken by 
the structure of the vacuum), a massless Goldstone boson is present, the PQ axion. Additionally, 
this model contains two CP-even Higgs bosons with masses of the order of the electroweak scale 
or lower. As in the MSSM, the heavy fields, one CP-odd, one CP-even and two charged Higgs 
bosons, are nearly degenerate. Fig.©. However, unlike the MSSM, the vacuum structure of 
the model constrains these heavy states to lie close in mass to ^ tan [3. Non-observation of the 
PQ axion rules out most of the allowed parameter space, only allowing scenarios with very low 
values of A and therefore very high VEVs for the new singlet field, which thereby do not provide 
a natural solution of the //-problem. 

If the PQ symmetry is slightly broken, the qualitative pattern for the particle spectrum re- 
mains intact, except that the lightest CP-odd state acquires a mass of the order of the electroweak 
scale. Thus the model contains a heavy CP-even and CP-odd Higgs boson plus two heavy charged 
Higgs bosons which are all nearly degenerate in mass, similarly to the MSSM. However, besides 
the standard light CP-even Higgs boson, the spectrum is complimented by a pair of CP-even 
and CP-odd Higgs bosons with masses below the electroweak symmetry breaking scale, Fig.(|H)). 
The LEP limits on the lightest Higgs boson mass allow some of the parameter space to be ruled 
out. However, since the couplings to the Z boson can be very much reduced, the NMSSM with 
a slightly broken PQ symmetry constitutes a valid scenario. Observing three light Higgs bosons, 
but no charged light Higgs boson, at future colliders would present an opportunity to distinguish 
the NMSSM with a softly broken PQ symmetry from the MSSM, even if the heavy states are in- 
accessible. It is therefore important that future colliders search for light scalar and pseudoscalar 
Higgs bosons with reduced couplings. 

In contrast, a strongly broken PQ symmetry, though disfavoured by the flow of the couplings 
from the GUT scale down to the electroweak scale, could provide extra moderately heavy Higgs 
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bosons, Fig.HlOj). which are only weakly coupled to the Z boson. Such decoupled scenarios would 
be more difficult to distinguish from the MSSM. 
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Appendix: An approximate solution 

The CP-even mass matrix of Eans. H23H28|) does not lend itself easily to obtaining analytic ex- 
pressions for the physical Higgs masses. However, reasonably simple expressions can be found by 
making an approximation, providing expressions which may be used to shed some light on the 
behaviour of the Higgs masses as the other parameters are varied. 

To construct this approximate solution in the scalar sector we regard both 1 / tan (3 and 1 /M^ 
as small parameters of magnitude e [Ma gauged by the generic electroweak scale]. Then, as 
long as neither k, A nor the other scales become too large, we observe a hierarchical structure in 
the CP-even mass matrix of the form: 




(67) 



where i? is a 2 x 2 matrix, C is a column vector and j4 is a scalar, all of order unity. 
Performing an auxiliary unitary transformation defined by the matrix. 




(68) 



with r = C/A, the mass matrix takes block diagonal form: 




(69) 
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where Mf^ are the entries of the CP-even mass matrix, Eans. ipSH^ -^^. Note that M^gZ-^l is 
actually 0{e^), so many terms can be neglected (and have been in the above). 

VM'^V^ is now easily diagonalized to give the approximate CP-even Higgs boson masses: 



M: 



H3 



M 



H. 



2/1 



MI2 + MI3 



Ml 



Ml + 



Mf, 
Ml, 



+ 4 



MI3 



Ml 




The diagonalization matrix is given by 

^10 \ 

R= cos Oh — sin 9h 
y sin 6h cos J 



with tan 0// 



(72) 



Combining R with the first rotation matrix V gives us the matrix O linking the states Si to the 
mass eigenstates (see Ean. (|lT|) : H = OS): 



O 



( 



\-My2Ml 
[Mfg sin Qh + Ml cos Qh\ I Ml 
\mI sin Oh - M^ cos Oh] /Ml 



Ml/M^ 



11 

cos 0H 

■ sin 9h 



\ 



Ml/Ml 
sm0H [1 - Mf^/lMl 
cos9h [l-Mf^/2Ml] J 



(73) 



This approximation works remarkably well, as can be seen in Fig.p2|). where both the exact 
(numerical) solution and the approximate solution, Eans. H7UH71|l . for the one-loop CP-even Higgs 
masses are shown as a function of M^i for the favoured parameter set of the slightly broken PQ 
symmetry, A = 0.3, k = 0.1, Vs = 3v, tan (3 = 3 and = —100 GeV, introduced in Sec. (3. 2). 
The exact and approximate solutions agree well as long as Ma does not become too small. 



^^In this section, we are concerned only with the CP-even squared mass matrix and the notation A! = M+ is 
understood. 
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Figure 12: The exact (numerical) solution and the approximate solution for the one-loop CP-even 
Higgs masses as a function of Ma for A = 0.3, k = 0.1, Vg = 3v, tan /? = 3 and = —100 Ge V. 
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